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' Abstract. We investigate whether it is possible to dynamically generate from 

classical transport theory the observed surprising R out » R S id e in Au+Au at 
y/s = 130 A GcV at RHIC [13, ||. We obtained covariant solutions to the 
Boltzmann transport equation via the MPC technique [ [|, for a wide range 
of partonic initial conditions and opacities. We demonstrate that there exist 
transport solutions that yield a freezeout distribution with R ou t < R S ide for 
K±_ >~ 1.5 GeV. These solutions correspond to continuous evaporation- like 
<^> freezeout, where the emission duration is comparable to the source size. Naively 

this would mean R ou t > R S ide- Nevertheless, our sources exhibit R ou t < R S ide 
because they are narrower in the 'out' than in the 'side' direction and, in 
addition, a positive x ou t — t correlation develops reducing R ou t further. 
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^ , 1. Introduction 

53 ■ 

Long-standing theoretical interests of the RMKI KFKI research group have been 
hydrodynamical, decoupling and transport phenomenaf |], [|, ||, Q). A recent theo- 
retical puzzle is why the pion source extracted via HBT interferometry in Au+Au 
at y/s — 130A GeV at RHIC has approximately identical 'out' and 'side' radii for 
average pair momenta K± < 0.7 GeV[[l], [^]. Most naively, this indicates very short 
pion emission duration, corresponding to almost instantaneous decoupling. 

The observed ratio R ou t / Rside ~ 0.9 — 1.1 contradicts theoretical expectations 
based on cylindrically symmetric and longitudinally boost invariant ideal hydrody- 
namics, indicating R ou t/ Rside significantly larger than one [ ^|, On the other 
hand, a hydrodynamical approach to the decoupling problem is inherently ques- 
tionable because in the decoupling region the mean free path exceeds the side of 
the system, while hydrodynamics is only valid if the mean free path is zero. As 
a workaround, hydrodynamics is usually supplemented with an arbitrary freeze- 
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out prescription postulating sudden freezeout on a three-dimensional hypersurface, 
which is typically parametrized by a constant freezeout temperature or energy den- 
sity. The unsatisfactory nature of this procedure has been realized recently and 
several improvements have been proposed [ ^| . 

A dynamical framework with a self-consistent freezeout is provided by classical 
transport theory [ 10, [H], O, |^, where the rate of interactions is given by the 



product of the density and the cross section A -1 = na. As the system expands, 
it rarifi.es, and the rate of particle interactions drops accordingly. Because the 
theory is not limited to equilibrium dynamics, it also allows for systematic studies 
of dissipative effects. Dissipation is crucial in order to reproduce[ |l| the observed 
saturation of elliptic flow at p± > 2 GeV[ [l4| |. 

The equilibrium assumption was partly relaxed in an earlier study [|15|, where 
hydrodynamical evolution was only assumed until hadronization at which point 
particles were fed into a hadronic cascade. The calculation gave R ou t / ' Rside ~ 1.2 — 
1.4, still above one. On the other hand, it is not clear whether the hadronic densities 
were large enough to support equilibrium on the freezeout hypersurface, which 
would be necessary to justify the assumed equilibrium evolution until the end of 
hadronization. In fact, antiprotons are known to scatter only once on average in this 
approach, indicating the need for a nonequilibrium approach before hadronization. 

In this study we investigate how nonequilibrium evolution in the partonic phase 
influences the pion distributions at freezeout. Utilizing the MPC numerical tech- 
nique [ |l6), we compute the covariant freezeout distributions for a variety of initial 
conditions for Au+Au at y/s — 130A GeV based on perturbative QCD, gluon satu- 
ration models, or a combination of both. Varying the transport opacity, which is an 
effective measure of the degree of collectivity [ we study under what conditions 
Rout < Rside is possible. 



2. Two-particle HBT interferometry 

Hanbury-Brown and Twiss interferometry is a truly remarkable technique to deduce 
space-time geometrical information purely via measuring momentum correlations. 
We review here only a few essential features and refer the reader to fiefs. [ [l?], [?], [ll| 
and references therein. 

Detailed studies of pion production by an ensemble of classical currents [ [l~7| 
showed that for a chaotic boson source p(x,t), in the limit of a large particle 
numbers, and also neglecting final state interactions, the two-particle correlation 
function is given by the space-time Fourier transform of the source as 



d x p(x, t)e l 



(1) 



Here = Pi — p^ is momentum difference, while = (pj* +p^)/2 is the average 
momentum of the boson pair. 

Unfortunately, Eq. (^) cannot be uniquely inverted to obtain the source dis- 
tribution, because the phase information is lost and, in addition, only a three- 
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dimensional projection of the four-dimensional correlation function is measurable 
because is constrained [|Tg|] by 

q^Kft = (ml — m 2 )/2 — (for identical particles) . (2) 

Despite the inversion problem, HBT measurements provide strong constraints on the 
freezeout distributions of particles and therefore on the possible dynamical scenarios 
in heavy-ion collisions. 

Conventionally, C(q,K) is expressed in terms of the 'out-side-long' variables, 
qo, <Zs, and q^. These are defined in the longitudinally boosted LCMS reference 
frame, where K z — with z or 'long' being the beam direction, 'out' the direction of 
K±, and 'side' the one orthogonal to both. In this frame, K£ CMS = (K°, K±, 0, 0), 
Qlcms = (qo,qs,QL), Xlcms = (t,x ,x s ,x L ) and Eq. © gives 

q^Xfj = -qo(xo - tK±/K°) - qs%s - Qlxl ■ (3) 

Clearly, the measurable projected correlation function is completely blind to the de- 
pendence of the source function on the fourth orthogonal coordinate t + xoK±/K°. 

Experimentally, the measured correlation function (after correcting for final 
state interactions) is fitted with a Gaussian, which for central collisions (b = 0) and 
at midrapidity is constrained by symmetry to be of the form 

C(q,K) = X(K)exp[-qlR 2 (K) - q 2 s R 2 s (K) - q 2 L R 2 L (K)\ . (4) 

Here Ro, Rs and Rl are the 'out', 'side', and 'long' HBT radii. For a perfectly 
Gaussian source, the correlation function is Gaussian and 

R 2 (K) = (Ax 2 ) K +v 2 ± (AP) K -2v ± (Ax At) K 

R 2 S (K) = (Ax 2 ) Kl and R 2 L (K) = (Ax 2 L ) K , (5) 

where v± = K±/ K a . Thus Rs and Rl have simple geometric interpretation as the 
'side' and 'long' widths of the distribution function, while Ro is a mixture of the 
'out' width, time spread, and the xo — t correlation. 

3. Covariant parton transport theory 

3.1. Transport equation 

We consider here, as in Refs. [ [ll], [l2[ |), the simplest but nonlinear form 
of Lorentz-covariant Boltzmann transport theory in which the on-shell phase space 
density f(x,p), evolves with an elastic 2^2 rate as 



pffy/i = S(x,fr) + JJJihU - hf2)W 12 ^6 4 ( Pl + P2 -P3-P4) ■ (6) 

234 
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Here W is the square of the scattering matrix element, the integrals are shorthands 
for J = J 7§^pj§7) where g is the number of internal degrees of freedom, while 

i 

fj = f(x,pj). The initial conditions are specified by the source function S(x,p). For 
our applications below, we neglect quark degrees of freedom and interpret f(x,p) as 
describing an ultrarelativistic massless gluon gas with g = 16 (8 colors, 2 helicities). 

The elastic gluon scattering matrix elements in dense parton systems were mod- 
eled with the isotropic form da e i/dt — ao(s)/s. The simplified angular dependence 
is justified by our previous study[ ||, where we showed that the covariant transport 
solutions do not depend explicitly on the differential cross section but only on the 
transport opacity 



0~tr 

X = W ~ Otr 



dzp (x +zh,T = ~J) . (7) 



o- e i J v c 

Here 

a tr (s) = J da el sin 2 9 cm = J dt ~^-^ ^1 - (8) 

is the transport cross section (in our case, o~tr = 2cro/3), while (n) is the average 
number of scatterings a parton undergoes. We neglected the weak logarithmic 
energy dependence of the total cross section and took a constant <To for simplicity. 
For a fixed nuclear geometry, a given transport opacity \ represents a whole class 
of initial conditions and partonic matrix elements, which is demonstrated by the 
approximate proportionality [ |j| x °c o~trdN g (To)/dr]. 

We solved Eq. (@) numerically via the MPC parton cascade algorithm [ |l6| |, 
which utilizes the parton subdivision technique [ EoL [ll| to obtain the correct, covari- 
ant solutions to the transport equation. Parton subdivision is essential to eliminate 
numerical artifacts caused by acausal (superluminal) propagation due to action at 
a distance [ H. For initial partonic densities expected at RHIC, the severe violation 
of Lorentz covariance in the naive cascade algorithm that employs no subdivision 
artificially reduces elliptic flow and heats up the p± spectra[ |[ 



3.2. Initial conditions 

We modeled central Au+Au collisions at RHIC with three different classes of initial 
conditions, corresponding to different assumptions on the initial transverse mo- 
mentum distribution and density distribution of partons. In all three cases, the 
evolution started from a longitudinally boost invariant Bjorken tube at proper time 
tq = 0.1 fm/c, with isotropic momentum distribution in the local rest frame and 
with uniform pseudorapidity i] = 1/2 log[(t+z)/(t—z)] distribution between |ry| < 5. 

The first class was essentially the same as the minijet initial conditions in Ref. 
[ [|. The initial transverse density distribution was proportional to the binary 
collision distribution for two Woods-Saxon distributions, while the p± distribution 
was a fit to the minijet distribution predicted by HIJING[ [l9| (without shadowing 
and jet quenching). To include the transport effect of the soft partonic component 
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present in HIJING in addition to the minijets, in the present study we increased 
the parton density by a factor five, to yield dN g (To)/drj = 1000. 

The second class was based on gluon saturation models. Here we assumed a 
uniform Bjorken cylinder with a radius i?o = 6 fm, dN g /drj = 1000 and a constant 
d 2 N/p±dp± distribution that vanishes for p± > Q sa t = 1.2 GeV. 

The third class was a combination of the first two. We took the minijet part 
with the original HIJING normalization (dN/drj — 210) and added four times as 
many saturated soft gluons. 

4. Results 

Using MPC, we computed the freezeout distribution numerically for central Au+Au 
at 130A GeV as a function of the average pair momentum and transport opacity, 
for all three classes of initial conditions. 

The freezeout distribution d 4 N/d 4 x was defined as the distribution of space- 
time coordinates for the point of last interaction of the test particles, with the 
strong assumption that this point is not affected by hadronization. We furthermore 
neglected resonance contributions to the pion yield. The same two hadronization 
schemes, parton-hadron duality and independent fragmentation, were applied as in 
Rcf. [|. 

Of our primary interest were the integrated distributions d 2 N/rdrdt, the equiva- 
lent of the usual hydrodynamical freezeout curves in the r—r plane, and d 2 N/dxodxs, 
which is the source projected onto the 'out-side' plane. We worked in the approx- 
imation px t ± w p2,_L ~ K±, which is valid up to corrections of 0(|(f|/|i?|). To 
get a simple estimate of the radii, we used Eq. (^|) even though our sources are 
non-Gaussian. A more detailed analysis of the radii based on the momentum-space 
correlation function is in progress. 

Figure [l] shows the freezeout distribution in r — t coordinates (r = \J Xq + x 2 s ) 
as a function of the pair momentum and transport opacity for the minijet initial 
condition with hadronization via parton-hadron duality. Unlike the sharp freezeout 
imposed in hydrodynamical models, the transport theory freezeout is a continuous[ 
|l2| ], evaporation-like process. For a given (nonzero) opacity, the larger the p± of 
the particle, the earlier it decouples and the closer it is to the surface of the nuclei. 
Thus, low-p± partons freeze out in the center at late times, while high-pj_ ones 
escape from the surface at early times. Furthermore, a larger opacity increases 
decoupling times, especially for low-pj_ particles. We observed the same features 
for the other two classes of initial conditions. 

Figure || shows the shape of the source in xo — x$ coordinates as a function 
of the pair momentum and transport opacity for the minijet initial condition with 
hadronization via parton-hadron duality. Results for the saturation and the 'com- 
bined' initial condition were very similar. The black crosses are positioned at the 
first moment of the distributions ((xo), (%s)), while their extension corresponds to 
the square root of the second moments (\J (Axq), (Air|)). As the opacity in- 
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creases, the source shrinks in the 'out' direction and the higher the p±, the more 
pronounced the reduction of (Ax ) is. From Eq. it is clear that this will help 
reduce Ro, at least if the other two terms do not change. On the other hand, we 
find that the source size in the 'side' direction is approximately independent of both 
the opacity and p±. Thus, surprisingly, R$ is unaffected by the strong collective 
dynamics. 
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Fig. 1. Freezeout r — t distribution as a fun<ji|gn] of transport opacity and px for 
the minijet initial condition with hadronization via parton-hadron duality. 



The mean (xo) does not coincide with the maximum of the distributions show- 
ing that the distributions are asymmetric in the 'out' direction, i.e., they are clearly 
non-Gaussian. In fact, the change in the shape of contour lines with increasing opac- 
ity is much more pronounced than the slow decrease of the second moment (Ax ) 
because there is a long tail in the negative 'out' direction. We expect that Ro 
extracted via a least-squared Gaussian fit would turn out to be somewhat smaller 



The decoupling problem at RHIC 



7 



than from Eq. (||). 
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Fig. 2. Freezeout xo — xs distribution aS(ji[fiinji c tion of transport opacity and K± 
for the minijet initial condition with hadronization via parton-hadron duality. 

The left column in Fig. |3| shows the dependence of the Ro/Rs and Rl/Ro 
ratios on transport opacity and K ^ for the minijet initial condition. For X = 0, both 
Ro and Rs are determined by the nuclear size to be as 3 fm, therefore Ro/Rs — 1- 
However, as the opacity increases, the ratio increases above one for K± < K Cl 
while decreases below one for Kj_ > K c , where the location of the turning point 
K c w 1.3 — 2 GeV depends on the hadronization scheme. 



The naive R Q 



R s + At expectation does not hold now because in the 



LCMS frame the source is not cylindrical (see Fig. |2j) and there are dynamical 
correlations. Figure || illustrates how the different terms in Eq. (|B|) contribute to 
Ro as a function of K±. For all three classes of initial conditions, the time spread 
(At 2 ) is comparable to the spatial width of the source ~ 3 fm, which has the effect 
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to increase Ro- However, this effect is compensated by the iarge positive xo — t 
correlation term, reducing Rq below R$- 
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Fig. 3. HBT radii ^igjj^ column) and ratios of radii (left columnX^s a function of 
transport opacity and K± for the minijet initial condition, with hadronization via 
parton-hadron duality (top and middle left) or independent fragmentation (bottom 
left and whole right column). 



The right column in Fig. £5] shows a comparison to the HBT radii measured at 
RHIC. In the transverse opacity range \ ~ — 8 we studied, the transport theory 
results are smaller than the observed Ro and R^. This is in sharp contrast to ideal 
hydrodynamics, which overpredicts both radii[ EI. The monotonic dependence of 
Rq and Rl on opacity indicates that a much better agreement with the data is 
possible at a larger but finite opacity x >^ 10. 

In Fig. ^, the reason for the anomalous small Rl for zero opacity is instan- 
taneous freezeout at the formation time t = tq. Since R\ w T 2 [A(n — y)] 2 , Rl 
depends on the decoupling time and the strength of the rj — y correlation. For our 
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thermally correlated initial condition [A (77 — y)] 2 ~ T/m± and To = 0.1 fm/c, which 
for x — gives Rl <~ 0.5 fm for all K± bins we studied. 

However, as the transport opacity increases, Rl becomes much larger because 
the decoupling time increases as is evident from Fig. [l]. The largest increase t/tq ~ 
R/tq ~ 50 is for low-p± partons, which freeze out latest. Thus the observed Rl(K±) 
carries important constraints on the dynamics of freezeout. For example, perfect 
inside-outside correlation, i.e., 77 = y as assumed in classical Yang-Mills approaches 
cannot be reconciled with the RHIC data because it gives Rl = 0. 

Finally, the biggest puzzle in Fig. ^ is why Rg(K±) ~ const « 3 fm contrary 
to the data and furthermore independent of the transport opacity. This suggests 
that Rg is independent of the degree of collective dynamics, i.e., it is determined 
solely by the initial nuclear geometry. The same problem with R$ being too small 
and insensitive to the dynamical assumptions has been observed in hydrodynamical 
calculations as well[ |). 




Fig. 4. HBTjjadjiiioj Rs an d magnitude of^ach term in Eq. (||) contributing to 
Ro as a function of K± for all three classes of initial conditions with hadronization 
via parton-hadron duality. 
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